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We consider a half-filled generalized Hubbard chain with electron-hole symmetric correlated hopping 
and on-site and nearest-neighbor repulsions U and V respectively. In addition to the insulating 
charge- and spin-density wave phases for large V and U respectively, we identify a phase with 
dominant superconducting correlations at large distances for small U and V . Using two Berry 
phases (one associated to the charge and the other to the spin degrees of freedom) as discrete order 
parameters, we construct a phase diagram for the three thermodynamic phases. 
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After the discovery of high- Tc systems, there has been 
an increase in the interest on purely electronic models 
exhibiting superconductivity, particularly those with cor- 
related hopping ||l|-[lOt. A hopping dependent on the oc- 
cupation of the two sites involved is to be expected in 
nature, and together with the on-site {U) and nearest- 
neighbor iy) repulsions, have been discussed already 
by Hubbard in 1963 [pdUj^ The resulting extended 
model (and a related one | [13| ) has been studied by sev- 
eral groups p^JI^ p7[| , with particular emphasis in phase 
transitions like the Mott metal-insulator one. At half- 
filling and when the hopping term which changes the 
double occupancy Iab = 0, Strack and VoUhardt ob- 
tained two regions of parameters for which the exact 
ground state in arbitrary dimension is a charge-density 
wave (CDW) or a spin-density wave (SDW) with max- 



imum order parameter |12|. These regions were later 
extended and a phase diagram constructed [ p^ . 

In one dimension (ID) the result is represented by the 
dashed line in Fig. 1 (b),(c). However, the nature of 
the third thermodynamic phase present in the diagram 
has not been clearly elucidated, no exact result exist for 
tAB 7^ 7^ and numerical studies of finite systems are 
not able to identify sharply the phase boundaries ||l^,|l8| . 
This difficulty is mainly due to the fact that all corre- 
lation functions of finite systems vary smoothly at the 
transition. 

On the other hand, the interest on applications of 
geometrical Berry phases has increased in recent years 
[^,|9-24|. In particular, Resta and Sorella have dis- 
covered that a Berry phase 7c (defined below) can be 
used as an order parameter for the ferroelectric transi- 
tion in a diatomic chain j2^. In ID problems with in- 
version symmetry (like ours) 7c can only attain the val- 
ues or TT (modulo 2ii) |jl^. This quantization of 7c 
in only two possible values does not allow to distinguish 
between three thermodynamic phases. Nevertheless, we 



define another Berry phase 7s, with similar features as 
7c, and show that there is a change in at least one of 
the topological numbers 7c/7i', Is/t^ at the transition be- 
tween any two of the thermodynamic phases. Therefore, 
each thermodynamic phase is univocally determined by 
a two-component Berry phase vector 7 = (7c, 7s). This 
fact allows us to define precisely all phase boundaries 
for any chain, small to obtain its ground state by the 
Lanczos method , and the phase diagram is obtained 
by finite-size scaling. In ID the nature of a thermody- 
namic phase is determined by the correlation functions 
(charge-charge, spin-spin or pair-pair) with the lowest 
decay at large distances. Calculation of the correlation 
exponents and other arguments allow to identify the ther- 
modynamic phase other than the CDW or SDW ones, 
as "superconducting" (S), in the sense that the system 
behaves as a Tomonaga-Luttinger liquid in which the su- 
perconducting pair-pair correlations are the greatest at 
large distances. 

The Hamiltonian for a ring of L sites with electron cre- 
ation operators satisfying arbitrary boundary conditions 
for both spins cj^^^. = e**''cjg. can be written in terms of 

the operators c]^ — e^^^"^" / ^c^-^ satisfying cj_|_^^ — 
as: 



JO- 

+tBBni^sni+i,s) + h.c.} 

+U 5Z "jT^ji + ^ H rij+i^rij^, (1) 

3 J0-0-' 

In 2D, this Hamiltonian has been derived as an effec- 
tive one-band model for cuprate superconductors [ p5| . 
We restrict the parameters to tAA = iBs = l,0< 
t-AB ^ 1, interpolating between the usual extended Hub- 
bard model Iab = 1 M, and the case tAB = 0, for which 
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some exact results are known ||,^,|l2|,|l5|,|6| , particularly 
the boundaries of the CDW and SDW phases and 
the exact ground state for V — Unfortunately 
this state is rather pathological, highly degenerate, and 
at half- filling it is insulating even for small U, in spite of 
the absence of a charge gap for U < A, for geometrical 
reasons [H,^- Thus, it is necessary to include the effect 
of tAB to draw general conclusions on the nature of the 
Mott or superconductor-insulator transition. 

Due to the periodic boundary conditions used, the 
many-body eigenstates of H for any fluxes $0-, can be 
classified according to their total wave vector K with 
KL/2iT integer (instead, the wave vectors K in the repre- 
sentation of the cl are shifted and cover all real values as 



the <i>cr are varied [||j2^j2|]). If an eigenstate |eif ($1, $x)} 
of i/(<I>|,$|) with a given K does not have an acciden- 
tal degeneracy with another eigenstate with the same K 
in the interval < $| = $x = $ < 27r, [|6| then this 
state can be followed adiabatically as $ is varied, and a 
many-body generalization of Zak's Berry phase [p^ can 
be defined for it |Q. For the explicit calculation, it is 
more convenient to use the numerically gauge invariant 
expression |2^,|2j], discretizing the interval < $ < 27r 
into TV points = 2TTi/N: 

IcicK) - - Jim {Im[ln(n^o2^eK($.„$.)|eK($.+i,$,+i)) 



(eK(<i>iV-i,<i>JV-i)|e^(27r)))]}, 



(2) 



where |e^(27r)) represents |ex(27r, 27r)) obtained directly 
from [6^(0,0)) [0: 



(27r)) = exp[z^^jc]^c,,]|eK(0,0)). 



(3) 



In practice iV = 10 points are more than enough to de- 
fine if 7c(e/f) is or ±7r with six digits accuracy. In 
few cases in which this was not the case, we increased 
TV until the same level of accuracy was obtained. Note 
that in contrast to previous cases ||2l| , p^ , |e/f ($,$)) is 
not necessarily the ground state. In particular, for most 
of the parameters inside the S phase, there is a crossing 
of the energy levels between K = and K = tt. As a 
consequence, the ground-state energy Eg{^, $) has local 
minima around $ = and $ = tt (see Fig. 2) and its 
form is suggestive of anomalous flux quantization (AFQ), 
i.e. Egijr) = Eg{0) in the thermodynamic limit, as ex- 
pected if superconducting correlations dominate at large 
distances |[||. 

In analogy to the fact that from the ground-state 
energy iJg($|,<i>x) one can calculate a charge stiff- 
ness (Drude weight) and a spin stiffness as Dc = 
{1/2)8"^ Eg{^,^)/d^^, Ds = (L/2)92£:g($,_$)/a$2^ 
for <I> = 0, respectively we define the spin Berry 
phase as that captured varying $| from to 27r, keeping 



$1 in the cycle: 



7s(e/f) 



— lim jlmfln 

{^f=o^{eK{<^^, -$,)|eK($.+i, -<i>.+i)) 
x(eK(<i>Ar-i,-<i>Ar-i)|e|f(27r)))]}, (4) 



27r 



e|f(2^)) = exp[i—^j(7c]^cj,] 16^(0,0)), 



(5) 
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where a — 1(-1) for spin up (down). While in general Tc 
is a measure of the polarization in the system [^-pT 24 
(a change of 7c in tt corresponds to the change of po- 
larization when one unit charge is displaced one lattice 
parameter), 7^ is a measure of the difference between the 
polarization of electrons with spin up and down. 

At this point it is interesting to note the effect of the 
unitary transformation. 



(6) 



on H. In the resulting Hamiltonian H' , <I>| and U change 
sign, while the term in V transforms to an Ising inter- 
action V' = V{nj^ — nji){nj+i^ — nj+ii). Assuming in 
the following L even, and calling |e^) the state obtained 
from |eif ) through the unitary transformation one has: 



-lc{eK) 



(7) 



The term tt is due to the fact that the transformation 
interchanges the vacuum with the fully polarized ferro- 
magnetic state for which "ic{F) — js{F) = tt. This 
can be easily seen from Eqs. (2)-(5), since all scalar prod- 
ucts inside the logarithm are 1, except the last one, which 
is -1. 

We define the Berry phases of a thermodynamic phase, 
as the Berry phases of the corresponding ground state. 
For those cases in which two energy levels with differ- 
ent quantum numbers as functions of $1, $x cross in the 
ground state, we have chosen the state of lowest mini- 
mum of the energy i?($| In all cases we have found 
that this minimum corresponds to X = $1 = $x = tt 
for L multiple of four, and to X = $x = $x = if i/4 
is a half integer. The chosen state has also the lowest 
energy averaged over the fluxes. The values of the Berry 
phases at the CDW and SDW thermodynamic phases 
are easy to determine using continuity arguments (7s , 7c 
cannot jump unless there is a phase transition). In the 
limit of very large U, for $| = <I>| = <I>, _ffis ecxuiv- 
alent to a Heisenberg Hamiltonian with exchange con- 
stant J = At\g/{U - V) independent of $. Then, all 
factors but the last in Eq. (2) are 1, the exponential in 
Eq. (3) introduces a factor -1, and 7c(SDW) — tt. A 
Hartree-Fock wave function leads to the same result, in 
agreement with previous findings for a diatomic system 
P4| . When $1 = —^i, the exchange terms in the large 
U limit depend on <I>|. However we do not expect a tran- 
sition if the Ising term (of the form of V') is increased 



2 



(artificially). This favors a Neel configuration and Eqs. 
(4) and (5) lead to 7s(SDW) = tt. Defining a vector 
7 — {ids), we can write 7(SDW) = (tt, tt). Similarly 
in the large V limit, the ground state is the CDW with 
maximum order parameter and from Eqs. (2)-(5), it eas- 
ily follows that 7(CDW) = (0,0). This result is consis- 
tent with Eq. (7) and with the changes in polarization 
for both spins when in the Neel SDW state, all electrons 
with spin up are displaced one lattice parameter to form a 
CDW state. These simple arguments cannot be extended 
to the S phase near U^V ^ 0. However for U = V = 
and ^ Iab 1, the ground state is non-degenerate 
for any fluxes $|, and therefore, it should be mapped 
onto itself by the transformation Eq. (6). Then, Eq. 
(7) implies that 7s(S) = tt + 7c(S) (modulo tt). Thus, 
before actually computing the Berry phases one knows 
that the 7 are different for the three phases. We obtain 
7(8) — (0,7r) numerically. In the non-interacting case 
U = F = 0, tAB = 1, 7 is ill defined due to a degeneracy 
of the ground state for some values of <i>| , $ (correspond- 
ing to a local maximum in the ground-state energy as a 
function of flux) and the same value of K . This is consis- 
tent with previous numerical calculations, which found 
that 7c in a metallic normal phase took arbitrary values 
and did not converge with the number of splittings N 
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FIG. 1. Phase diagram in the U, V plane for different values 
of tAB- Solid (empty) circles denote the points at which 7c 
(7s) jumps in tt. Dashed line corresponds to the exact result 
for tAB — + [^. Dashed-dotted line in (a) corresponds to 
V = U/2. 

The phase diagram for different values of tAB is shown 
in Fig. I. Because of the unitary transformation Eq. 
(6), the phase diagram is also that of H' interchang- 
ing SDW and CDW phases. The results were obtained 
calculating the boundaries for all chains of even length 
L < 10 and extrapolating the results using the proce- 
dure followed in Ref. For tAB = 1, the S phase 
is absent and the CDW-SDW boundary is very similar 
to that found previously from the structure factors ob- 
tained with Monte Carlo simulations |l^. However, this 
method has an error ^ 0.05 in the critical value of V. 
There is in addition a systematic error due to the dis- 
cretization of imaginary times which is difficult to esti- 
mate. In our case, from the difference in the extrapolated 
values using up to 8 sites, or up to 10 sites, we estimate 
the error to be one order of magnitude smaller. As long 
as tAB < 1, a superconducting phase appears, near the 
CDW-SDW boundary for i^s = 1, if f/ < 4, particu- 
larly for small U and V. This can be understood in part 
by the fact that the energy of CDW and SDW states 
increase when tAB decreases. This fact is particularly 
clear for large V or large U for which the ground-state 
energy per site is Ecdw = U/2 — 2t\g/{3V — U) and 
EsDW = V — itAB In 2/([/ — F) respectively, and agrees 
roughly with the result of a generalized Hartree-Fock de- 
coupling of the two- and three-body terms present in H 
0. The exact resuhs for tAB = @,|,|l|,^|6) , show 
that for small tAB, U and V , it is energetically convenient 
to destroy the CDW or SDW phase to take advantage of 
the kinetic energy terms in tAA and tBB- However, the 
fact that the S phase is more stable near the line in which 
the CDW and SDW are degenerate is noticeable and un- 
expected in these previous approaches. This is probably 
related to the fact that the transition between the CDW 
and SDW phases for tAB = 1 is second order for U <Uc 
and first-order for U > Uc, where the position of the 
tricritical point Uc ^ 3 ± 1 [|l8|. The S phase seems to 
develop around the line of the second-order phase tran- 
sition as long as tAB < 1- Unfortunately, our method 
can not determine the order of the transition and a more 
precise value of Uc can not be given. 

As tAB decreases further, the region occupied by the S 
phase increases and tends to the exact result for tAB 
|l6| . For small tAB the finite-size effects become more im- 
portant for the S-SDW phase boundary. For V — and 
tAB = 0.1, the extrapolation using chains up to 10 sites 
gives Us-SDW = 3.901, while extrapolation including 
L — 12 gives a value smaller by 0.016. 
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FIG. 2. Ground-state energy as a function of the flux 
= = ^ for U = 1.5, V = 0, tAB = 0.5 and several 
values of L. 

Previous numerical studies for finite i^B, well inside 
the S phase have shown the presence of mobile carriers 
, however the attempts to determine if the S phase 
is superconducting or not gave contradictory results 
On one hand, it has been established that for tAB < 1, 
and small U and V, the system behaves as a Tomonaga- 
Luttinger liquid with dominant superconducting correla- 
tions at large distances. On the other hand, no signs 
of AFQ were found for C/ = ^ = 0, and a general- 
ized Hartree-Fock Bardeen-Cooper-Schrieffer (BCS) de- 
coupling of the two- and three-body terms, shows that 
an ordinary singlet BCS superconducting solution is un- 
stable. However, we find that this decoupling, for small 
U, V, leads naturally to a p- wave triplet-pairing supercon- 
ducting ground state. A similar BCS solution exists in 
2D, which should be valid for small U, V and 1 — tAB- In 
addition, since the dependence of the ground-state energy 
with the flux is calculated simultaneously with the Berry 
phase 7c, we found that in most of the S phase there are 
signs of AFQ. An example is given in Fig. 2. The fact 
that no signals of AFQ are found ior U = V = is prob- 
ably a finite-size effect. Thus, all calculations are con- 
sistent with the superconducting nature of the S phase. 



Furthermore, for V = 0, the correlation exponent Kp > 1 
for U < Us-SDW, while for U > Us-sdw a gap opens 
in the charge sector Calculations of Kp for Iab = 0.6 
in rings of 8, 10 and 12 sites gave Us-sdw = 2.05 ±0.05, 
in excellent agreement with the result 2.11 obtained with 
the Berry phases (see Fig. 1 (b)). The triplet nature of 
the S phase is supported by our numerical calculations 
of pair-pair correlation functions for L = 10 and the fact 
that for a singlet S phase one expects 7 = (0, 0). 

In summary, we have shown that in one dimension at 
half filling, the extended Hubbard model with correlated 
hopping H (Eq. (1)), and the model H' obtained from 
H through the unitary transformation Eq. (6), present 
a phase "S" with dominating superconducting correla- 
tions at large distances for < tAB < tAA ~ tsB 
and small U and V. This thermodynamic phase, as 
well as the insulating CDW and SDW phases , which 
are the ground state for large V and large U respec- 
tively, are characterized by different quantized values of 
7 — (7c, 7s), where 7c, 7s are Berry phases associated 
with charge and spin degrees of freedom respectively: 
7(S) = (0,7r), 7(CDW) = (0,0), 7(SDW) = (7r,7r). The 
finite jumps in 7 at the phase boundaries allowed us to 
define them sharply for each finite chain, and to obtain 
an accurate phase diagram through finite-size scaling. 
The use of topological quantum numbers to character- 
ize the different thermodynamic phases should be useful 
in general for non-trivial highly correlated models, for 
which correlation functions in finite systems do not show 
sharp transitions and exact solutions in the thermody- 
namic limit are very rare. 
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